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ABSTRACT 

y  /> 

For  the  exponential  life  distribution  model  and 
any  prior  distribution  for  the  failure  rate 
parameter,  the  predictive  distribution  has  a  de¬ 
creasing  failure  rate.  'Ue-give^^a  Bayesian  ex¬ 
planation  of  why  this  is  logically  reasonable. 


A  BAYESIAN  EXPLANATI(»I  OF  AN  APPARENT  FAILURE  RATE  PARADOX 


by 

Richard  E.  Barlow 

For  many  devices,  wear-out  does  not  seem  likely — at  least  not  within 
the  time  frame  of  Interest.  For  such  devices,  an  exponential  life  dis¬ 
tribution  model  may  be  used;  l.e.,  for  random  lifetime  X  , 

P[X  >  X  I  X]  -  e“^*  . 

However,  Information  about  X  must  come  from  data  and/or  engineering 
Judgment.  In  the  exponential  case  and  relative  to  X  ,  the  data  can  be 
summarized  by  k  ,  the  observed  number  of  failures  and  T  ,  the  observed 
total  time  on  test.  Given  (k,T)  ,  suppose  a  posterior  density  Is  cal¬ 
culated  based  on  a  prior  ir(X)  .  The  predictive  survival  probability  to 
age  X  Is  then 


F(x  I  k,T)  -  I  P[X  >  X  I  X]tt(X  1  k,T)dX  .  (1) 

*1) 

It  Is  well  known  that  for  all  priors,  ir  ,  the  predictive  distribution  has 
a  decreasing  failure  rate  function  In  x  [cf.  Theorem  4.7,  Page  103, 

Barlow  and  Proschan  (1981)]  when  it  exists.  However,  at  first  glance,  this 
result  seems  highly  unreasonable.  If  the  device  does  not  wear  out,  why 
should  we  predict  Its  future  life  by  a  model  [the  predictive  distribution] 
which  actually  has  a  decreasing  failure  rate? 

Denote  this  predictive  failure  rate  function  by  r(x  |  k,T)  where. 


using  (1), 


2 


r(x  I  k,T) 


Def 


I 


Xe"^\(X  I  k,T)dX 


e"^*TT{X  1  k,T)dX 


(2) 


Figure  1  Illustrates  our  apparent  paradox. 


FIGURE  1 

COMPARISON  OF  PREDICTIVE  FAILURE 
RATE  AND  MODEL  FAILURE  RATE,  X 


A  Bayesian  explanation  of  Figure  1  may  be  instructive.  Note  that  from  (2), 


r(0  I  k,T) 


Xtt(X  I  k,T)dX 


E[X  I  k,T] 


(3) 


so  the  posterior  mean  for  X  estimates  r(0  |  k,T)  .  Now  suppose  a  new 
device  (exchangeable  with  our  sample  devices)  were  to  survive  x  hours, 
then  from  (2) 


'I 
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3 


r(x 


! 


I  k,T)dX 


k.T)  - 


OB 

I 


e~^\(X  I  k,T)dX 


But,  by  Bayes  theorem 


-Xx 


Tr(X  I  k,T,X  >  x)  -  — — -  I 


/ 


e~^%(X  I  k,T)dX 


so  that 


r(x  1  k.T)  -  EtX  I  k,T,X  >  x) 
“  E[X  1  k.T  +  x]  . 


(4) 


Hence  our  estimate  for  r(x  |  k.T)  based  on  our  updated  information,  namely 
that  another  device  has  survived  time  x  ,  produces  an  estimate  for  X 
different  from  r(0  |  k.T)  . 

It  Is  proved  In  Barlow  and  Proschan  (1979)  that  for  g  nondecreasing, 

CD 

J*  g(X)iT(X  I  k,T)dX  Is  decreasing  in  T  .  Expression  (3)  Is  the  special 
0 

case  g(X)  -  X  so  that  E[X  |  k.T]  ^  ElX  |  k.T  +  x]  for  all  priors. 


Example: 


baxa-l--bX 


If  it(X  I  a.b)  ■  - -  ,  the  natural  conjugate  prior,  then 

I  ■  (b  S 


and 


J 


4 


r(x  I  k,T) 


a 

b  +  T  +  X 


which  Is  obviously  decreasing  In  x  . 


Conclusion; 

The  source  of  the  confusion  concerns  the  difference  between  a  model 
based  on  prior  assumptions,  such  as  constant  failure  rate,  and  a  ^predictive 
distribution  based  on  current  Information.  We  would  not  use  the  predic¬ 
tive  distribution  (1)  as  our  model,  since  we  believe  In  constant  failure 
rate.  However,  were  we  asked  to  predict  whether  or  not  a  new  device 
(exchangeable  with  the  sample  devices)  would  survive  to  age  x  ,  we  would 
use  (1)  based  on  current  Information. 

On  the  other  hand,  were  we  to  characterize  a  new  set  of  say,  m  , 
devices,  exchangeable  with  our  sample,  we  would  say  that  each  has  an  ex¬ 
ponential  life  distribution  with  expected  failure  rate  E[X  |  k,T]  and 
irCX  I  k,T)  would  fully  measure  our  uncertainty  about  that  failure  rate. 
However,  were  we  asked  to  predict  how  xoany  will  fall  In  time  period  [0,t]  , 
we  would  use  the  predictive  distribution  and  calculate 

mEj^[l  -  1  k,T]  -  mF(t  1  k,T)  , 
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